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Method and device for estimating carrier frequency offset of 

subscriber terminal 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention relates to Code Division Multiple Access (CDMA) 
for mobile communications, more particularly to a method and device for 
estimating carrier frequency offset (CFO) in subscriber terminals in Time 
Division Synchronous Code Division Multiple Access (TD-SCDMA) 
system. 

Description of the Related Art 

In the TD-SCDMA system, both base stations and terminals are 
supposed to send and receive data with nominal carrier frequency. However, 
there is always some deviation between the actual transmitting and receiving 
frequency and the nominal frequency because of the technology limitation of 
the devices. In a base station, the frequency accuracy of an oscillator, in 
which a carrier frequency offset in a base station required by a TD-SCDMA 
system is less than 0.05PPM, may meet to the standard requirement, since 
limitations subjected by a base station from temperature, size, power 
consumption, cost and the like are smaller. In a subscriber terminal, however, < ' 
a carrier frequency offset in a subscriber terminal required by the 
TD-SCDMA system is less than 0.1 PPM, the frequency accuracy of a 
selected oscillator may not meet the standard requirement, due to the 
limitations subjected by the subscriber terminal from temperature, size, 
power consumption, cost and the like are greater. Consequently, carrier 
frequency offsets need to be estimated in a subscriber terminal by using a 
known synchronous downlink pilot signal (SYNC_DL) sent by base stations, 



and according to the estimated value of carrier frequency offsets to perform 
automatic frequency control (AFC), so that carrier frequency offsets in the 
subscriber terminal may meet the standard system requirements. 

In a method for estimating carrier frequency offset in prior art, 
SYNC_DL signals only send by a single base station are utilized in a 
subscriber terminal to estimate the carrier frequency offset. The base station 
has one or more antennas to transmit signals with the same carrier frequency. 
FIG. 1 is a schematic diagram showing a carrier frequency offset estimation 
system having a plurality of base stations in the prior art. 

In the TD-SCDMA system, a dedicated downlink pilot time slot 
(DwPTS) is set, a base station sends SYNC_DL signals at DwPTS of each 
frame. Codes by which a SYNC_DL signal is consisted of are used to 
distinguish adjacent cells using the same carrier frequency, that is, the 
adjacent cells with the same carrier frequency are distinguished by the 
SYNCDL codes. There are 32 different SYNCDL codes defined in the 
standard, length of each set of SYNC_DL codes is 64 chips. During the 
initial search of the cells, the subscriber terminal search the SYNCDL codes 
at a specific frequency point used in a cell where it currently expects to 
resided in to implement the downlink synchronization, and then estimate the 
carrier frequency offset based on the SYNC JDL codes and adjust the carrier 
frequency according to the estimated carrier frequency offset value to meet 
the system requirements. 

Provided that the system requires for carrier frequency offset estimation 
and rough adjusting, adjusted carrier frequency offset in a subscriber 
terminal is less than 1 kHz, then the subscriber terminal has only received, at 
a corresponding frequency point, the SYNC_DL signals transmitted by one 
base station or received S YNC__DL signals of the cell currently expected to 
resided in are much greater than the SYNC_DL signals of other cells using a 
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same frequency. If the signal power ratio is over 6dB, it is possible for the 
adjusted carrier frequency offset of a subscriber terminal to meet to the 
system requirement. However, if the power difference between the 
simultaneously received SYNC_DL signals respectively in two cells with the 
same frequency is small, the SYNCJDL signals transmitted by one cell will 
interfere those transmitted by another cell. For example, when the power 
ratios of multiple S YNC JDL signals received by a subscriber terminal at the 
edge of several cells are within 3dB, there is about 10-20% probability with 
the existing methods for carrier frequency offset estimation and adjustment 
that the carrier frequency offset of the subscriber terminal can't meet the 
requirement of the system. 

In the TD-SCDMA system, because of the short length of a SYNCDL 
signal, which is only 64 chips (50(is), the carrier frequency offset after 
adjustment permitted by the standard is still greater, for example, less than 
1kHz. Even so, as shown in FIG 1, when the subscriber terminal 
simultaneously receives SYNC_DL signals from several base stations, and 
whose power differences are small such as within 3dB by using an existing 
method for carrier frequency offset estimation, which only uses the received 
SYNC_DL signal in a cell (base station) where it currently expects to reside 
in, i.e., the highest power SYNC_DL signal, to make a carrier frequency 
offset estimation, there is about 10-20% probability that the carrier frequency 
offset in the subscriber terminal after automatic frequency control can't be 
below 1kHz. In this case, the initiate search for a cell by a subscriber 
terminal will not be implemented and probability of success for the initiate 
search for a cell will be affected. 

Summary of the Invention 

In view of this, an object of the present invention is to provide a method 
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and a device for estimating carrier frequency offset in a subscriber terminal 
to make the carrier frequency offset in a subscriber terminal to meet the 
system requirements and thus improve the probability of success for the 
initiate search for a cell. 

To achieve the above object, the present invention provides a method 
for estimating carrier frequency offset in a subscriber terminal, said method 
comprises steps: 

A. determining number of effective base stations from which more than 
one signals are received by a subscriber terminal and main path 
positions of each signal; 

B. combining the signals of each station corresponding to said number 
of effective base stations based on the main path positions obtained in step 
A; 

C. calculating a rough estimation value of the carrier frequency offset 
based on combined signal in step B. 

Said determining number of effective base stations from which more 
than one signals are received by a subscriber terminal in step A comprises 
steps: 

Al. calculating peak power value of each signal received by a 
subscriber terminal, and selecting the peak power values of predefined 
number of base stations from higher to lower; 

A2. determining the number of effective base stations from predefined 
signals are received by the subscriber terminal by comparing the ratio of the 
highest peak power value from the order in step Al to the subsequent peak 
power values with the given threshold. 

Said signals are synchronous downlink pilot signals, and said step Al 
further comprises steps: 

All. shift multiple correlating a local synchronous downlink pilot code 
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and a received synchronous downlink pilot signal results in a power value of 
the synchronous downlink pilot signals received by the subscriber terminal ; 

A 12. determining peak power values corresponding to each of the 
synchronous downlink pilot codes. 

Said step A2 further comprises steps : 

A21. numbering the peak power values ordered from the highest to the 
lowest and setting a current sequence number as predefined number of the 

base stations ; 

A22. determining whether the highest peak power value and a peak 
power value corresponding to the current sequence number are greater than 
the given threshold, if so, setting the number of effective base stations from 
which the signals are received by a subscriber terminal as the value of the 
current sequence number, otherwise, the current sequence number decreases 
by one and returns back to step A22. 

The method further comprises a step before said step A : reading vector 

data of 128 chips while receiving synchronous downlink pilot signals at the 
beginning of a downlink pilot time slot. 

The method further comprises a step before said step B: multi-path 
combining signals of each base station. 

Said step of multi-path combining signals of each base station 

comprises steps : 

beginning from a point of previously predetermined number of the peak 
power value, reading data of synchronous downlink pilot signals at a point 
which is 2 times of the predetermined value added length of said 
synchronous downlink pilot code; 

performing Max Ratio Combination after eliminating phase difference 
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between symbols of multi-path synchronous downlink pilot signal with 
different time delay and the phase difference of delay path. 

Said step B of combining the signals of each station corresponding to 
the number of effective base stations is: equal gain combining or weight 
combining signals of each base station corresponding to said effective 
base station number to obtain an combined signal sequence. 

Said step C is to obtain a rough estimating value of the carrier frequency 
offset according to the phase difference between two symbols spaced by a 
defined distant in said combined signal sequence. 

Said step C further comprises: estimating carrier frequency offset for a 
predefined times, and then averaging them to get a carrier frequency offset 
estimation. 

Said step C is to sum up the phase differences between two symbols 
spaced by a defined distant in said combined signal sequence, and then 
computing the phase angle to get the carrier frequency offset estimation. 

The invention also provides a device for estimating carrier frequency 
offset, which comprises at least: 

a decision module for determining effective base station number from 
which signals are received by a subscriber terminal and a main path 
position of signal transmitted from each base station based on the signals 
received by a subscriber terminal, and then outputting the number of the 
effective base station and the main path position of each signal to an 

combining module ; 

an combining module for combining the signals from each base station 
corresponding to the number of effective base stations based on the main 
path position of signals and then outputting the combined signals to a 
carrier frequency offset acquiring module; 

a carrier frequency offset acquiring module for obtaining a rough 
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estimating value of the carrier frequency offset based on the combined 
signals. 

Said device further comprises a multi-path combining module for 
multi-path combining the signals of each base station, and then outputting 
the multi-path combined signal to the combining module, if the base station 
number is greater than 1 

According to a solution of the invention, when a subscriber terminal 
simultaneously receives SYNC_DL signals sent by a plurality of base 
stations (more than one) and their power differences are small, it is possible 
to make the carrier frequency offset in a subscriber terminal, after automatic 
frequency control, to meet the system requirement and thus improve the 
probability of success for the initiate search for a cell. 

The method for estimating carrier frequency offset provided by the 
invention is especially suitable to a condition that the strengths and SIR of 
SYNC_DL signals received by a mobile user from several base stations are 
approximate among them, i.e., a mobile user is located at the joint edge of 
several cells, and with the method for carrier frequency offset estimation 
provided by the invention, the probability of success for the initiate search 
for a cell can be improved greatly. Even when the SYNC_DL signal received 
by a mobile user from one base station is much stronger and those received 
from other base stations are not so strong, the accuracy of adjustment for 
carrier frequency offset also can be improved by using the method for 
estimating carrier frequency offset provided by the invention. With respect to 
a method utilizing the SYNC_DL signals from one base station, the closer of 
the SYNCDL signal strength and SIR from several base stations, the more 
improvement the accuracy of carrier frequency offset adjustment. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



FIG. 1 shows a schematic diagram of a prior art system having several 
base stations for estimating carrier frequency offset ; 

FIG. 2 shows a flow chart for estimating carrier frequency offset in a 
subscriber terminal according to the invention ; 

FIG 3 is a detailed flow chart showing the determination of the 
SYNC_DL signals simultaneously received by a subscriber terminal from 
several base stations with the same frequency and the arriving time from its 
main path in the steps as shown in FIG. 2; 

FIG 4 is a flow chart principally showing the multi-path combination of 

the SYNC_DL signals received from one base station ; 

FIGS shows the simulation outcome according to an embodiment of the 
invention. 

DETAILED DESCRIPTION 

Hereinafter, detailed description will be made to a method and device 
for estimating carrier frequency offset in a subscriber terminal in 
TD-SCDMA system according to the embodiments of the invention, with 
reference to accompanying drawings, so that the invention will be well 
understood by those skilled in the art. 

FIG 2 shows a flow chart for estimating carrier frequency offset in a 
subscriber terminal according to the invention, and FIG. 3 is a detailed flow 
chart showing SYNC_DL signals simultaneously received by a subscriber 
terminal from several effective base stations with the same frequency and 
arriving time from respective main path determined in step 201 shown in FIG. 
2. 
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The first step, in step 201, determine the SYNC_DL signals 
simultaneously received by a subscriber terminal from several effective base 
stations with the same frequency and the respective arriving time of their 
main path. 

FIG. 3 illustrates the detailed flow chart in step 201. In step 301, after 
the synchronization of DwPTS in a subscriber terminal, 128 chips of data r 
are read beginning with the start of DwPTS, where the kth element is 
denoted as r k , k-0, 1, 127. 

Then, in step 302, shift complex correlation with the data r received at 
DwPTS to 32 groups of local SYNC_DL codes respectively, shifting from 0 
to 63, i.e., calculate the correlation power values of each point. 

In step 303, the correlated power values for K frames are averaged 
(K>1). 

The signal powers for two successive frames arrived at a subscriber 
terminal may be different because of fading, etc,. To ensure the accuracy of 
the data obtained in a subscriber terminal, it is preferable to calculate the 
correlation power by using more frames of data and then average them. 

Then, in step 304, averaged correlation power values corresponding to 
SYNC_DL codes of each group in 32 groups are fined out; the arriving time 
for the SYNCDL signals corresponding to the SYNC_DL codes are 
determined by the SYNC_DL codes, i.e., the arriving time for the main path 
of the SYNC_DL signals corresponding to the SYNCJDL codes; and, shift 
complex correlation with data r to the SYNC_DL codes to obtain the peak 
value of correlation power, which indicates the synchronization of received 
signals and the corresponding SYNCJDL codes. 

In step 305, the peak values of the averaged correlation power are 
ordered in a falling sequence, and the first N max ones are selected, denoted 

as p max( n ) i n== l> 2, N max . N max is a predefined value, and represents the 
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max number of base stations in which a subscriber terminal receives the 
SYNCDL signals from different base stations simultaneously. 

In step 306, since a subscriber terminal is not able to know how many 
SYNCDL signals are received simultaneously from base stations when it 
receives SYNC_DL signals, thus it is needed to determine the number N of 
effective base stations from which the S YNC DL signals have been received 
simultaneously, based on a threshold, set the threshold V T , and firstly decide 
whether the N max th peak value of the averaged correlation power meets 
following condition: 
P (1) 

maxW_ > y 

P (N ) 

max V max / 

If yes, N = N max ; otherwise, continue to decide whether the (N max -l)th 
peak value of the averaged correlation power values meets following 
condition: 

SnaxO) > y 

P (N -1) 

max V- 1 * max / 

If yes, N = - 1 ; otherwise, the rest may be deduced by analogy, 

until the number N of effective base stations, from which a subscriber 
terminal receives the SYNCJDL signals simultaneously, is determined, 
where l<N<N max . The SYNC_DL codes sequence number and the peak 

values point corresponding to each of the N base stations are recorded as 
Sync _ID n and Peak n respectively , where n= 1,2, ...,N. 

In the embodiment of the invention, N max =3, and Vt=4. 
The second step, in step 202, combine multipath signals of SYNCJDL 
received from each effective base station respectively. 

Then, the complex SYNC_DL codes numbered as Sync_ID n are 

operated, and the complex values of SYNC_DL codes are marked as s n , 
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where the kth element is s n k , k=0, 1, . .., 63. 64+2M points of data are read 
beginning from M chips ahead of Peak n in the DwPTS of the received 
signals, marked as a vector e , where the kth element is e k , k=l, 2, 
63+2M. 

e is shift complex correlated with s n , and the result cor n is obtained : 

63 

cor n(i =X e i + k < k . i = 0 5 1, ... ? 2M , 

k=0 

then the SYNC_DL signals of each paths are combined and a vector 
U n is obtained : 

2M 

U n , k =JXk -s^ -corV , k=0 3 1, 63, 

where, s* nk indicates a conjugate of s nk , and cor n . indicates a 
conjugate of cor n . . 

The physical meaning of the formula is to perform Max Ratio 
Combining of multi-path SYNC_DL signals with different time delay after 
eliminating the phase differences both between the symbols and between the 
time delay paths, in which, received signals e multiplied by the conjugate 
of s n corresponds to eliminating the phase difference between symbols, 

and multiplied by the conjugate of cor n corresponds to eliminating the 

phase difference between time delay paths. 

The process of step 202 of multi-path combining SYNC_DL signals 
received from a subscriber terminal will be further described with reference 
to FIG 4, as shown in FIG.4. 

In step 401, 64+2M points of data are read beginning from M chips 
ahead of Peak n in the SYNC_DL signals received at DwPTS, i.e., 64+2M 
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received data is read beginning from M chips ahead of the position of the 
main path, and a complex sequence of 64+2M is obtained, e(k) , k=0, 1, . . 

63+2M , marked as a vector e , where 64 is the length of the SYNC_DL 
code. 

In step 402, the SYNC_DL codes numbered SYNCJD are complex 
operated and the complex SYNC_DL codes are obtained , i.e., the local 

complex SYNC_DL code sequence s(k) , k=0, 1, 63 , marked as s n 0 

Then, in step 406, conjugating operation of each element of the 
sequence s(k) is performed, and a sequence s*(k) is obtained. 

Because e(k+l) includes a phase shift, in order to eliminate the phase 

shift between subsequent two phases, in step 403, letd(0,k) = e(k) ,k=0, 1, . .., 

63 ; in step 404 , let d(l,k) = e(k + l) , k=0, 1, ...,63 , and in step 405 , 

let d(2M,k) = e(k + 2M) , k=0, 1, ...,63. 

In step 407 , each element of the sequence d(i,k) ( i = 0,1,- --,2M , k=0, 

1, 63 ) obtained in steps 403-405 is multiplied by the corresponding 
element of the sequence s*(k) obtained in step 406 respectively. 

In step 408 , all elements of the sequence obtained in step 407 are 

summed up and then conjugated. In steps 403 -step 408, shift multiplex 
correlation of e and s* are performed , and a result cor n is obtained as 

follows : 

63 

cor„,i =2Xk<k , i=0 5 1, ...,2M , 
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here, e i+k is the d(i,k) and the e(i+k). 

In step 409 , the result of shift multiplex correlation is multiplied by all 

elements of the sequence, then, in step 410, all 2M+1 sequences, i.e., 
multi-path combining the SYNC_DL signal of each path, are summed up and 

a vector U n is calculated as follows : 

2M 

U n , k =£ e i+k < k corV , k=0, 1, 63 , 

i=0 

The value of M is integrally dependent on both the time delay extension 
of the received signals and the effect of noise. In the embodiment of the 
invention, M = 4. 

The third step, in step 203, multi-path combined SYNC_DL signals of 
each base station are further combined. 

combining the multi-path combined S YNC DL signals of each effective 
base station is to combine the calculation results of the SYNC_DL codes, 

and equal gain combination is used : 

U k =IXic , k=0,l,...,63 ; 

n=i 

and a weighting combination may also be used : 
U k =2X k .P n , k=0, 1 63 , 

n=l 

where the weighting coefficients P n can be obtained by the following 
formula : 

2M 2 

P„=S| cor ".i| • i = 0,l,-,2M , n = l,2 > -,N , 

i=0 

Because there is coherent property among the SYNC_DL signals from 
each effective base station while noncoherent property among noises, it is 
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possible to improve the SNR (or SIR) of available signals by combining the 
multi-path combined SYNC_DL signals from each effective base station so 
that the SYNCDL signals received by a subscriber terminal may be 
efficiently utilized. 

Based on theory formula, the sequence combined by the SYNCDL 
signals from each effective base station can be approximately represented as 

follows : 

U k =A k -exp(j-2^-Af -k-T c +^ 0 )TfN k , k=0, 1, 63 , 

where, U k is the kth symbol value in the combined SYNC_DL signal 

sequence, k=0, 1, 63 ; 

A k is magnitude of combined SYNCDL signals valid in the kth 

symbol ; 

exp(x) indicates exponent operation of complex x; 

Af is carrier frequency offset, which is the exact value in theory ; 

T c is pulse width of each chip; 

<z) 0 is initial phase of the SYNC__DL signals received by a subscriber 
terminal, range of which is between O-ln ; 

N k is the complex noise value or interference noise value in the kth 
symbol o 

The fourth step, in step 204, rough estimation values of carrier 
frequency offset are calculated by using the combined signals. 

Based on the theoretical formula in step 203, it can be seen from the 
expression of U k that in the combined SYNCJDL signal sequence, a 
variation rate of phase between symbols is associated with a carrier 
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frequency offset Af . Therefore, a carrier frequency offset can be roughly 
estimated by using the combined SYNCJ3L signal sequence. In one 
embodiment of the invention, 

I T ^Zangle(U;-U k+L ) 
2;rLT c PS 

Where Af is a rough estimation of the carrier frequency offset, T c is 
an pulse width of each chip, and angle (x) is an operation of the phase angle 

of complex x. Values P and L will meet formula: P+L<64 , and P>1 , L>1. 

Above formula indicates the effect of noise (or interference) can be 
minimized by using a phase difference between two symbols spaced by L to 
estimate the carrier frequency offset, and then averaging the estimated values 
for P times . 

In another embodiment of the invention, 

Af= 1 angle[X(U;-U k , L )] 
In L- l c k=0 

Where Af is rough estimated value of the carrier frequency offset, T c 
is pulse width of each chip, and angle (x) presents the operation of the phase 

angle of complex x. Values P and L will meet formula: P+L<64 , and P>1 , 
L>1. 

In order to minimize calculating amount of repeatedly calculating phase 
angle, the above formula indicates to firstly sum up the phase differences 
between two symbols spaced L, and then calculate phase angles to estimate 
the carrier frequency offset. 

Due to the phases of two consecutive symbols are close together, an 
angle jitter of the multiplied phase angle will be greater. Therefore, it is 
preferable to multiply two symbols spaced L to avoid the angle jitter to be 
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more great. 

In the embodiment of the invention, P = L = 32 0 

The invention also provides an device for estimating carrier frequency 
offset in a subscriber terminal, which corresponding to the steps shown in 
FIG.2, comprising a decision module for determining number of base 
stations, from which the subscriber terminal receives SYNC_DL signal, 
based on the SYNCJDL signal received by the subscriber terminal and the 
main path positions of the SYNC_DL signals from each base station, and 
then outputting the number of effective base stations and the main path 
position of each signal to an combinting module; an combining module for 
multi-path combining the SYNCJ3L signal from each effective base station 
based on the main path position of a signal , and then combining the 
multi-path combined SYNC_DL signal of each effective base station. The 
number of SYNC_DL signal in each effective base station to be combined is 
the number of effective base stations determined in the decision module, and 
then outputting the combined SYNC_DL signal to a frequency offset 
acquiring module; and a frequency offset acquiring module for estimating 
carrier frequency offset of the subscriber terminal according to the combined 
SYNC signals received. The device may further comprises a multi-path 
combining module for multi-path combining the signals from each effective 
base station and outputting the multi-path combined signal to the combining 
module. 

FIG.5 shows the simulation outcome of the embodiment of the invention. 
In this simulation, we assume that the subscriber terminal (UE) receives 
simultaneously the SYNCJDL signal from 3 base stations having same 
frequency, and the relative power received from three base stations are OdB, 
-0.5dB and -ldB respectively, with corresponding SNR of 3dB, 2.5dB and 
2dB. The deviation of a carrier frequency from a nominal frequency of these 
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3 base stations are uniformly distributed within (-100Hz, 100Hz). The 
relative path time delays of the three SYNC signals received by the 
subscriber terminal are distributed uniformly within 0-4 chips, with time 
delay resolution of 1/8 chip. The initial carrier frequency offset of UE is 
distributed uniformly within (-10kHz, 10kHz). The channel condition is an 
additive white Gaussian noise (AWGN) channel. 

UE carrier frequency offset is roughly adjusted by using the method 
provided by the invention. Parameters are selected as follows: K=4, V T =4, 
M=4. The simulation outcome is shown in Figure 5. N rnax =l indicates that 
the carrier frequency offset is estimated by using the SYNC signal from one 
base station, N max =2 indicates that the carrier frequency offset is estimated 
by using the SYNC signal from two base stations, at most, having same 
frequency, and N max =3 indicates that the carrier frequency offset is estimated 
by using the SYNC signal from three base stations, at most, having a same 
frequency. It can be seen from FIGS, when UE simultaneously receives the 
SYNC signals from 3 base stations having same frequency and transmission 
powers from each base station is not much different from each other, if the 
carrier frequency offset is estimated by using the SYNC signal only from 
one base station (N max =l), then there is still about 27% of carrier frequency 
offset to be over 1kHz after AFC. But if the carrier frequency offset is 
estimated by using the SYNC signal from three base stations having same 
frequency, it is possible to keep the carrier frequency offset under 1kHz and 
the adjustment accuracy thereof will be increased greatly. 

In summery, the description above is only the preferred embodiments of 
the invention but will not limit the protection scope of the invention. 
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